Quantum Phase Transitions Without Thermodynamic Limits 
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A new microcanonical equilibrium state is introduced for quantum systems with finite-dimensional 
state spaces. Equilibrium is characterised by a uniform distribution on a level surface of the ex- 
pectation value of the Hamiltonian. The distinguishing feature of the proposed equilibrium state 
is that the corresponding density of states is a continuous function of the energy, and hence ther- 
modynamic functions are well defined for finite quantum systems. The density of states, however, 
is not in general an analytic function. It is demonstrated that generic quantum systems therefore 
exhibit second-order (continuous) phase transitions at finite temperatures. 
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The derivation of phase transitions in quantum sta- 
tistical mechanics typically requires the introduction of 
a thermodynamic limit, in which the number of degrees 
of freedom of the system approaches infinity. This limit 
is needed because the free energy of a finite system is 
analytic in the temperature. But phase transitions are 
associated with the breakdown of the analyticity of ther- 
modynamic functions such as the free energy. Hence in 
the canonical framework the thermodynamic limit is re- 
quired to generate phase transitions. Although the ex- 
istence of this limit has been shown for various systems 
(see, e.g., 0), the procedure can hardly be regarded as 
providing an adequate description of critical phenomena. 

One can consider, alternatively, a derivation based on 
the microcanonical ensemble. The usual construction 
of this ensemble @ is to define the entropy by setting 
S = ks In tie j where tie is the number of energy levels in 
a small interval [E,E + AE}. The temperature is then 
obtained from the thermodynamic relation TdS = dE. 
This approach, however, is not well formulated because 
(a) it relies on the introduction of an arbitrary energy 
band AE, and (b) the entropy is a discontinuous func- 
tion of the energy. To resolve these difficulties, a scheme 
for taking the thermodynamic limit in the microcanoni- 
cal framework was introduced in For finite systems, 
however, the difficulties have remained unresolved. 

The purpose of this paper is to demonstrate the fol- 
lowing: (i) if the microcanonical density of states is de- 
fined in terms of the relative volume, in the space of pure 
quantum states, occupied by the states associated with a 
given energy expectation E, then the entropy of a finite- 
dimensional quantum system is a continuous function of 
E, and the temperature of the system is well defined; 
and (ii) the density of states so obtained is in general not 
analytic, and thus for generic quantum systems predicts 
the existence of second-order phase transitions, without 
the consideration of thermodynamic limits. 

It is remarkable in this connection that similar types of 
second-order transitions have been observed recently for 



classical spin systems, for which the associated configura- 
tion space possesses a nontrivial topological structure 0] . 

The paper is organised as follows. We begin with the 
analysis of an idealised quantum gas to motivate the in- 
troduction of a new microcanonical distribution. This 
leads to a natural definition of the density of states fl(E). 
Unlike the number of microstates Ue, the microcanonical 
density Q(E) is continuous in E. As a consequence, we 
are able to determine the energy, temperature, and spe- 
cific heat of elementary quantum systems, and work out 
their properties. In particular, we demonstrate that in 
the case of an ideal gas of quantum particles, each particle 
being described by a finite-dimensional state-space, the 
system exhibits a second-order phase transition, where 
the specific heat decreases abruptly. 

Ideal gas model. Let us consider a system that con- 
sists of a large number N of identical quantum parti- 
cles (for simplicity we ignore issues associated with spin- 
statistics). We write -fftotai for the Hamiltonian of the 
composite system, and Hi (i — 1,2, ... , N) for the Hamil- 
tonians of the individual constituents of the system. The 
interactions between the constituents are assumed to 
be weak, and hence to a good approximation we have 
-Htotai = J2iLi Hi- We also assume that the constituents 
are approximately independent and thus disentangled, so 
that the wave function for the composite system is ap- 
proximated by a product state. 

If the system as a whole is in isolation, then for equi- 
librium we demand that the total energy of the com- 
posite system should be fixed at some value £?total- In 
other words, we have (iftotal) = -^totai- It follows that 
Yli=i{Hi) = -Etotai- Now consider the result of a hy- 
pothetical measurement of the energy of one of the con- 
stituents. In equilibrium, owing to the effects of the weak 
interactions, the state of each constituent should be such 
that, on average, the result of an energy measurement 
should be the same. That is to say, in equilibrium, the 
state of each constituent should be such that the expec- 
tation value of the energy is the same. Therefore, writing 
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E = N^ 1 .©total) we conclude that in equilibrium the gas 
has the property that (H) — E. That is to say, the 
state of each constituent must lie on the energy surface 
Ee in the pure-state manifold for that constituent. Since 
N is large, this will ensure that the uncertainty in the 
total energy of the composite system, as a fraction of 
the expectation of the total energy, is vanishingly small. 
Indeed, it follows from the Chebyshev inequality that 
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for any choice of a; > 0. Therefore, for large N the energy 
uncertainty of the composite system is negligible. 

For convenience, we can describe the distribution of 
the various constituent pure states, on their respective 
energy surfaces, as if we were considering a probability 
measure on the energy surface Ee of a single constituent. 
In reality, we have a large number of approximately in- 
dependent constituents; but owing to the fact that the 
respective state spaces are isomorphic we can represent 
the behaviour of the aggregate system with the specifi- 
cation of a probability distribution on the energy surface 
of a single "representative" constituent. 

Microcanonical equilibrium. In equilibrium, the 
distribution is uniform on the energy surface, since the 
equilibrium distribution should maximise an appropri- 
ate entropy functional on the set of possible probabil- 
ity distributions on £ E . From a physical point of view 
we can argue that the constituents of the gas approach 
an equilibrium as follows: On the one hand, weak ex- 
changes of energy result in all the states settling on or 
close to the energy surface; on the other hand, the inter- 
actions will induce an effectively random perturbation in 
the Schrodinger dynamics of each constituent, causing it 
to undergo a Brownian motion on Ee that in the long run 
induces uniformity in the distribution on E E - We con- 
clude that the equilibrium configuration of a quantum 
gas is represented by a uniform measure on an energy 
surface of a representative constituent of the gas. 

The theory of the quantum microcanonical equilibrium 
state presented here is analogous in many respects to 
the symplectic formulation of the classical microcanoni- 
cal ensemble described in There is, however, a subtle 
difference. Classically, the uncertainty in the energy is 
fully characterised by the statistical distribution over the 
phase space, and for a microcanonical distribution with 
support on a level surface of the Hamiltonian the en- 
ergy uncertainty vanishes. Quantum mechanically, how- 
ever, although the statistical contribution to the energy 
variance vanishes, there remains an additional purely 
quantum-mechanical contribution. Hence, although the 
energy uncertainty for the composite system is negligible 
for large N, the energy uncertainties of the constituents 
will not in general vanish. An expression for AH will be 
given in equation Q below. 



Density of states. To describe the equilibrium rep- 
resented by a uniform distribution on the energy surface 
Ee, it is convenient to use the symplectic formulation of 
quantum mechanics. Let TL denote the Hilbert space of 
states associated with a constituent. We assume that the 
dimension of TL is n + 1 . The space of rays through the 
origin of TL is a manifold T equipped with a metric and 
a symplectic structure. The expectation of the Hamilto- 
nian along a given ray of TL then defines a Hamiltonian 
function H(ip) = (i/j\Hi\i/j) / (ip\4>) ° n I\ where the ray 
if) e r corresponds to the equivalence class \ip) ~ A|t/>), 
A G C\0. The Schrodinger evolution on TL is a symplec- 
tic flow on r, and hence we may regard T as the quan- 
tum phase space. Our approach to quantum statistical 
mechanics thus unifies two independent lines of enquiry, 
each of which has attracted attention in recent years: the 
first of these is the "geometric" or "dynamical systems" 
approach to quantum mechanics, which takes the sym- 
plectic structure of the space of pure states as its start- 
ing point ; and the second of these is the probabilistic 
approach to the foundations of quantum statistical me- 
chanics in which the space of probability distributions on 
the space of pure states plays a primary role Q • 

The level surface £ e in V is defined by H(ip) = E. The 
entropy associated with the corresponding microcanoni- 
cal distribution is S(E) = fcBlnf2(B), where 



Sl(E) = J S(H(t(j) - E)dV r . 



(2) 



Here dVr denotes the volume element on T. In a mi- 
crocanonical equilibrium the temperature is determined 
intrinsically by the thermodynamic relation TdS = dE 1 
which implies that k B T = Vt(E)/9!(E), where tt'(E) = 
dfl(E)/dE. Since the density of states fl(E) is differen- 
tiable, the temperature is well-defined. Other thermo- 
dynamic quantities can likewise be precisely determined. 
For example, the specific heat C(T) = dE/dT is given 

b y c = k B {n') 2 /[{n') 2 -nn"}. 

Consider a large system composed of two independent 
parts, each in a state of equilibrium. Each subsystem 
is thus described by a microcanonical state with sup- 
port on the Segre variety corresponding to disentangled 
subsystem states. Let us write Qi(Ex) and ^2(^2) for 
the associated state densities, where E\ and E 2 are the 
initial energies of the two systems. Now imagine that 
the two systems interact weakly for a period of time, 
during which energy is exchanged, following which the 
systems become independent again, each in a state of 
equilibrium. As a consequence of the interaction the 
state densities of the systems will now be given by ex- 
pressions of the form Qi(Ei + e) and ^(JSa — e), for 
some value of the exchanged energy e. The value of e can 
be determined by the requirement that the total entropy 
S(E) = k B hx\p,\(E\+e)Q,i(E2— e)} should be maximised. 
A short calculation shows that this condition is satisfied 
if and only if e is such that the temperatures of the two 
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systems are equal. This argument shows that the defi- 
nition of temperature that we have chosen is a natural 
one, and is physically consistent with the principles of 
equilibrium thermodynamics. 

Phase transitions. The quantum microcanonical 
ensemble introduced here is applicable to any isolated 
finite-dimensional quantum system for which the ideal 
gas approximation is valid. The volume integral in PJl 
can be calculated by lifting the integration from T to H 
and imposing the constraint that the norm of \ip) is unity. 
Then we can write: 
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E)dV n , (3) 



where dV-n is the volume element of H. Making use of 
the standard Fourier integral representation for the delta 
function, and diagonalising the Hamiltonian, we find that 
© reduces to a series of Gaussian integrals (see [a| for 
details). Performing the ^-integration we then obtain the 
following integral representation for the density of states: 



Q(E) = (-mf 
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The two integrals appearing here correspond to the delta- 
functions associated with the energy constraint H(ip) = 
E and the norm constraint (i/j\i/j) = 1. Carrying out the 
integration we find that the density of states is given by 
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where 1{a} denotes the indicator function (1{a> = 1 if 
A is true, and otherwise). In (J3J) we let m denote the 
number of distinct eigenvalues Ej (j — 1,2, ... , m), and 
we let Sj denote the multiplicity associated with the en- 
ergy Ej. Thus Y^JLi $j = n>+ 1. In the nondegenerate 
case, for which Sj = 1 for j = 1, 2, 



(-7T)" 



n+1 



. , m 



n+1 



(n-1) 



we have 



L {E k >E} 

Ei — E/, 



(6) 



With these expressions at hand we proceed now to ex- 
amine some explicit examples. 

Nondegenerate spectra. In the case of a Hamilto- 
nian with a nondegenerate spectrum of the form E^ = 
e(k — 1), k = 1, 2, . . . , n + 1, where e is a fixed unit of 
energy, the density of states © reduces to 
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We see that Q,(E) is a polynomial of degree n — 1 in each 
interval E £ [Ej , Ej + i), and that for all values of E it 
is at least n — 2 times differentiable. In Fig. ^ we plot 
Q(E) for several values of n. For a system in equilibrium 
the accessible values of E are those for which Q'(E) > 0. 
States for which Sl'(E) < have "negative temperature" 
in the sense of Ramsey . 



n(E) 




FIG. 1: Density of states Q(E) for nondegenerate (n + 1)- 
level systems for n = 3, 4, 5, with energy eigenvalues Ek = k 
(k = 0, . . . , n). Q(E) is n — 2 times differentiable in E. 

The structure of the space of pure states in quan- 
tum mechanics is intricate, even for relatively elemen- 
tary systems. In particular, as the value of the energy 
changes, the topological structure of the e nerg y surface 
undergoes a transition at each eigenvalue For ex- 

ample, in the case of a nondegenerate three-level system, 
the topology of the energy surface changes according to: 
Point — » S 3 — » S 1 x — > S 3 — > Point, as the energy is 
raised from -E m ; n to E max (R^ denotes a two-plane com- 
pactified into S 2 at a point corresponding to the interme- 
diate eigenstate). These structural changes in the energy 
surfaces induce a corresponding nontrivial behaviour in 
the thermodynamic functions. 

As an illustration we consider a four-level system and 
compute the specific heat as a function of temperature. 
The result is shown in Fig. [2 where we observe that 
the specific heat drops abruptly from 2fcs to \ks at the 
critical temperature T c defined by ksT c = \e. Therefore, 
this system exhibits a second-order phase transition, in 
this case at the critical energy E c = e. This example 
shows that the relationships between phase transitions 
and topology discovered recently in classical statistical 
mechanics carry over to the quantum domain where, 
arguably, they may play an even more basic role. 

For a system with a larger number of nondegenerate 
eigenstates, the specific heat also increases abruptly as T 
is reduced. In this case the specific heat is continuous, 
and the discontinuity is in a higher-order derivative of the 
energy. For a system with n+1 nondegenerate energy 
eigenvalues, the (n — l)-th derivative of the energy with 
respect to the temperature has a discontinuity. The phe- 
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nomenon of a continuous phase transition is generic, and 
is also observed if the eigenvalue spacing is not uniform. 

Degenerate spectra. In a system with a degener- 
ate spectrum, the phase transition can be enhanced. In 
particular, the volume of Ee increases more rapidly as E 
approaches the first energy level from below, if this level 
is degenerate. This leads to a more abrupt drop in the 
specific heat (Fig.[3J). 



at E — E max we note that the first moment of £l(E) is 
given by the integral of H (x) over T. Hence by use of a 
trace identity obtained in [l2| we have 
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FIG. 2: Specific heat for a nondegenerate four-level system 
(dotted line, n = 3, Ej = 0, 1, 2, 3), a four-level system having 
a degenerate first excited state (dashed line, n = 4, Ej — 
0, 1, 1, 2, 3), and a quantum Ising chain (solid line, J = 1/4, 
B = 1). In the quantum Ising chain, we have C(T) ~ (T — 
T c )~ 2 away from T c , whereas in the vicinity of T c we have 
C(T) ~ (T - T c )~ 13 for T > T c . (We set k B = 1 here.) 

As an example, we consider a quantum Lenz-Ising fer- 
romagnetic chain with three spins. The Hamiltonian is 

H = -./ELi cr M +1 ~ -SELi ^ where a z denotes 
the third Pauli matrix for spin fc, and J, B are con- 
stants. We have in mind a gas of weakly interacting 
molecules, each modelled by a strongly-interacting quan- 
tum Ising chain. The eigenvalues of the Hamiltonian arc 
Ei = —3 J - 35, £ 2 ,3,4 = J - B, E 5>6 ,7 = J + B, and 
E$ = — 3J + 3B. As the temperature is reduced, the spe- 
cific heat grows rapidly in the vicinity of the critical point 
T c = (2 J + B) /iks (Fig. , where the system exhibits a 
discontinuity in the second derivative of the specific heat. 
We note that when B is small the critical temperature is 
close to that of the classical mean-field Ising model. 

Density matrix and energy uncertainty. Finally, 
we show the existence of a natural energy band associ- 
ated with the quantum microcanonical distribution. The 
microcanonical density matrix for the energy E is 

fiE = -^r J S(H(^)-E)il(fl)) dVr. (8) 

Here TL(ip) = (ip\il>) denotes the projection oper- 

ator onto the state S TL corresponding to the point 
■ijj G r. The squared energy uncertainty is (AH) 2 — 
tr(fiEH 2 ) — [tr(fiEH)] 2 . A calculation then shows that 

{AH)2 = ^ Se {R ~ U) Q{U) dU ' (9) 
where H = tr(H)/(n + 1) denotes the uniform average 
of the energy eigenvalues. To check that AH vanishes 



u Q(u) du 



H(ip) dVr = —r H. (10) 



However, the integral of Q,(E) is the volume 7r"/n! of T, 
and the desired result follows. Using the explicit formulae 
obtained earlier for il(E) we are then able to calculate the 
energy uncertainty associated with the equilibrium state 
of a finite quantum system. It remains to be seen whether 
the new ensemble can be put to the test in some definitive 
way, and in particular whether the phase transitions it 
predicts actually correspond to observable phenomena. 
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